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ABSTRACT 

Disk formation in magnetized cloud cores is hindered by magnetic braking. Previous work has 
shown tliat for reahstic levels of core magnetization, the magnetic field suppresses the formation of 
rotationally supported disks during the protostellar mass accretion phase of low-mass star formation 
both in the ideal MHD limit and in the presence of ambipolar diffusion for typical rates of cosmic ray 
ionization. Additional effects, such as ohmic dissipation, the Hall effect, and protostellar outflow, are 
needed to weaken the magnetic braking and enable the formation of persistent, rotationally supported, 
protostellar disks. In this paper, we first demonstrate that the classic microscopic resistivity is not 
large enough to enable disk formation by itself. We then experiment with a set of enhanced values 
for the resistivity in the range rj = 10^^-10^^ cm^ s~^. We find that a value of order 10^^ cm^ is 
needed to enable the formation of a 10^ AU-scale Keplerian disk; the value depends somewhat on the 
degree of core magnetization. The required resistivity is a few orders of magnitude larger than the 
classic microscopic values. Whether it can be achieved naturally during protostellar collapse remains 
to be determined. 

Subject headings: accretion, accretion disks — ISM: clouds — ISM: magnetic fields — magnetohydro- 
dynamics (MHD) 



1. INTRODUCTION 
1.1. Magnetic Braking and Disk Formation 

Circumstellar disks play a central role in both star and 
planet formation. The disks are thought to be the con- 
duit throug h which Sun-like stars assem ble the bulk of 
their mass ( Shu, Adams, fc Lizano|[T987 ) . They are also 
the birthplace of planets Despite ye ars of active research 
(see reviews by Bodenheimer 1995 for early work), we 
still do not have answers to some of the most basic ques- 
tions, such as "when and how does a disk appear in the 
process of star formation?" 

Observationally, it is difficult to separate out the pro- 
tostellar disk from the surrounding envelope during the 
early, embedded phase of star formati on with the cur- 
rent g eneration of interferometers (e.g., J0rgensen et al. 
20071, although the situation will improve " drastically in 



a tew years, when ALMA comes online (van Dishoeck 



& J0rgensen 


20081 


istence of a 


large 



20081. Enoch et al. (2009) interred the ex- 



(more than O.IMq) disk around the protostar Serpens 
FIRS 1, based mostly on CARMA dust continuum ob- 
servations. Without a detailed knowledge of the central 
stellar mass and gas kinematics on small (arcsec) scales, 
it is difficult to determine whether the inferred structure 
is a rotationally supported disk, a magnetically induced 
"pseudodisk," or simply the inner part of an envelope 
that is dynamically more complicated than the Ulrichjs 



( 1976 1 rotating, collapsing sphere adopted by [Enoch et 
aE ( |2U09 ) for the protostellar envelope. Theoretically 



Hlsk formation in dense cores of molecular clouds mag- 
netized to the observed level (corresponding to a mean 
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dimensionless mass-to-flux ratio A ~ a few; Troland & 
Crutcher|2008l ) is controlled by magnetic braking, which 
has been difficult to quantify until recently. 

Detailed calculations have now shown that the mag- 
netic braking is apparently so efficient as to inhibit 
the formation of rotationally supported disks during the 
main protostellar accretion phase for a moderate level 
of core magnetization, as long as the magne tic field is 
frozen in the matter (i.e., the ideal MHD limit; I Allen, Li 
& Shu 2003; GaUi et al. 2006; Price & Bate 2007; Hen^ 
nebelle fc Fromang,2008; .Mellon fc Li,2008; see, however, 
Hennebelle fc Ciardi||2009| for a difi'erent view, an d also 



Machida et al.|20 05land |Banerjee fc Pudritz|2006| for the 



evolution of rotating, magnetized cores prior to the main 
accretion phase). In order for a rotationally supported 
disk to exist around a rapidly accreting protostar, the 
magnetic braking must be weakened one way or another. 
The most obvious possibility is through non-ideal MHD 
effects, which are expected to be important in lightly 
ionized star-forming cores. However, we have recently 
shown that ambipolar diffusion, the best studied non- 
ideal MHD effect in star formation, does not sufficiently 
weaken the braking to allow rotationally supported disks 
to form for realistic levels of core magnetization and ion- 
ization; in some cases, the mag netic braki ng is even en- 



see 



han ced dMellon fc L i 2009, Krasn opolsky fc Konigl 2002 



(]2l 



Basu fc Mousch ovias 1995, Hos king"fc Whitwort 

|2004| and [Puffin fc Pudritz, 2009, who considered the ef- 
fects of ambipolar diffusion on magnetic braking and core 
fragmentation during the earlier phase of prestellar core 
evolution) . This motivates us to examine whether two of 
the remaining non-ideal MHD effects, ohmic dissipation 
and the Hall effect, can weaken the coupling between 
the magnetic field and the bulk neutral matter (and thus 
the strength of magnetic braking) enough to enable the 
formation of rotationally supported disks. This paper 
concentrates on the ohmic dissipation. 
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1.2. Electrical Resistivity of Cloud Cores 
The resistivity rj is related to the electrical conductivity 



three orders of magnitude, see Figure 1 of 



(T through 



a 

Ana 



(1) 



For lightly ionized molecular gas in dense cores, the elec- 
tric conductivity is usually dominated by the contribu- 
tion from electrons, although the grain contribution can 
dominate at high densities when large amounts of small 
grains are present (e.g., Wardlc & Ng 1999). To obtain a 
rough estimate for the conductivity, we will concentrate 
on the contribution from the electrons, which is given by 



mep{crv)e--ii2 



= 9.0 X lO^^Xe (s-1) , (2) 



where 



adopted a value (av) 



2.8TOHn-H2 E^nd Xe = rie/nfja- We have 



2 X 10-y cm'^ 



for a 10 K 

gas from Pinto & Galli ( 2008[). This va lue is somewhat 
smaller than tha t used in Uraine et ai. ( 1983 ) and War- 
dle fc Ng| ( [l999| , which is 2.6 x IQ-^^Js^ but some- 
what larger than that of Sano et al. ( 2000[ ), which is 
1.3 X 10~^cm^s~^. It increases slowly with tempera- 
ture, roughly as T^/^, yielding a lower resistivity at a 
higher temperature. Since our calculation is assumed to 
be isothermal, the temperature dependence is not used. 
Substituting equation ([2| into equation ([T]) yields 

7.9 X 102 _ 
rj = ( cm^ s ^) . (3) 

The fractional ionization Xf, can be computed from de- 
tailed chemical networks. At the relatively high densities 
that are relevant for disk formation, it depends on the 
grain s ize distribution. For the s tandard [Mathis, Rumpl, 
fc Nord sieck ( 1977[ ) distribution, Nakano, Nishi, fc Ume- 
bayashi| (|2002| computed the ionization traction over a 
large range of number density (their Figure 1). Over 
the critical range of density ^ 10^-10^^ cm""^ where the 
transition between the rapidly infalling pseudodisl|j^and 
a potential Keplerian disk is expected to occur, theirac- 
tion can be fitted roughly by 



10" 



^14 



10^° cm- 



which corresponds to a resistivity of 



V ■ 



1.7 X 10 



17 



10~i3gcm 



(4) 



(cm^s-^). (5) 



The fit of fractional ionization deviates substantially 
from the computed values at low densities. However, the 
deviation is expected to be inconsequential because the 
low density envelope is sufficiently ionized to be close to 
the ideal MHD limit (ignoring ambipolar diff usion, which 
is investigated separ ately elsewhere, see e.g., 
sky fc Konigll [20021 and 



to either equation 



2]ar 



2009 



Krasnopol- 



accordmg 



Mellon fc Li| 

or the computed values. If only 
relatively large grams of size (of order 0.1 /xm) exist, 
the ionization fraction would be higher (by as much as 

By "pseudodisk" we mean in this paper a structure formed 
by the magnetically induced flattening of infalling; material due to 
anisotropic magnetic forces, as in ^Galli fc Shu^^l993J . 



ashi 
ardle 



"fc Nakano 1990 or compare Figures 1 and 3 of 

fSi Ng||1999 ), making the resistivity smaller. Equation 

t should therefore be viewed as a rough indication of 
I plausible upper range of the "classic" (as opposed 
to "anomalous" ) resistivity, especially in view of the fact 
that the actual temperature in the high density region 
of disk formation is expected to be higher than the 10 K 
adopted in our calculation, which would lower the resis- 
tivity. It may b e enhanced, however, throu gh anomalous 



processes (e.g., Norman fc Heyvaerts| 1985 ) 

The rest of the paper is organized as follows. In §2, 
we describe the setup of the disk formation problem, fo- 
cusing on the protostellar mass accretion phase, when a 
central object has formed. For illustrative purposes, we 
will concentrate on a particular set of core parameters 
for which a 10^ AU-scale disk is formed at a representa- 
tive time t = 10^2 s (or about 3 x 10^ yr) in the absence 
of magnetic braking, and a magnetic field that is strong 
enough to suppress the formation of a rotationally sup- 
ported disk in the ideal MHD limit. We consider in §3 
whether the illustrative classic resistivity given in equa- 
tion ([5]) can weaken the magnetic braking enough to en- 
able disk formation, and find that the answer is no. In §4, 
we explore the question of how large the resistivity must 
be in order for the 10^ AU-scale rotationally supported 
disk to reappear in our model problem. It turns out to be 
of order 10^^ cm^ s""'^, a few orders of magnitude larger 
than the class ic value expected on the 10 AU-scale, as 
anticipated by Shu et al. ( 2006 1 based on the concept of 
the so-called Ohm sphere, although our value is one to 
two orders of magnitudes smaller than theirs, depending 
on the field strength. The numerical results and their 
astrophysical implications are discussed in the last sec- 
tion, §6. 

2. PROBLEM SETUP 

Our primary goal is to illustrate the effects of ohmic 
dissipation on magnetic braking and disk formation. For 
this purpose, it is desirable to have as small a numer- 
ical magnetic diffusivity as possible, which demands a 
high spatial resolution. High resolution, however, puts a 
stringent limit on the computation time step dt, which 
is proportional to the square of the (smallest) grid size 
for the explicit method that we use to treat the ohmic 
dissipation. Although subcycling can be used to allevi- 
ate the time step problem substantially, a large number 
of cycles are still needed to reach a reasonable protostel- 
lar accretion time, say, 10^^ s, making the treatment of 
self-gravity difficult. For this initial exploration of the 
effects of ohmic dissipation, we have decided to set up a 
simplified disk formation problem, with a rotating, non- 
self-gravitating envelope falling onto a central object of 
fixed mass Af* . It will turn out that a large "anomalous" 
resistivity is needed to save the disk, an elaborate disk 
formation model is probably not warranted at this early 
stage of the investigation. 

Specifically, we solve the equations of resistive MHD, 
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dp 
dt 
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dB 

~dt 



47rJ = cV X B 



V X (v X B - 77V X B) 



(8) 
(9) 



under the assumption of axisymmetry, and adopt a 
spherical polar coordinate system {r,9,(j)), with a stel- 
lar object of M* = O.5M0 at the origin such that the 
gravitational potential is $g = —GM^/r, neglecting the 
self-gravity of the gas. Self-gravity will be included in a 
follow-up work, where we will evolve the magnetized core 
from its formation to collapse in a self-consistent manner. 
For this initial study, we aim to capture the essence of 
the problem of magnetic suppression of disk formation 
in the simplest possible way, by adopting uniform dis- 
tributions for the initial gas density and magnetic field, 
and to study the effects of a new ingredient, resistivity, 
in this simplest problem. 

At time t = 0, we fill the computation domain be- 
tween Ti = 1.5 X 10^''' cm and Vo = 1.5 x 10^^ cm with 
a uniform density po — 1-4 x 10~"'^^gcm~^ (correspond- 
ing to ri(H2) — 3.0 X lO^cm"'^), so that the total enve- 
lope mass is IMq. For simplicity, we assume that the 
gas stays isothermal, with an isothermal sound speed 
^ = (p/pY^^ — 2 X lO^'cms^^ (corresponding to a tem- 
perature of about 10 K). For the initial rotation, we 
adopt the following prescription: 



(10) 



where we choose = 2 x 10^ cms ^, n7 is the cylindri- 
cal radius, and -cUc = 3 x 10^^ cm. The softening of the 
rotational profile inside the cylindrical radius Wc is to 
prevent the angular speed from becoming singular near 
the rotation axis. Since the inner radius of our com- 
putational domain is 10 AU, we are concerned with the 
formation or suppression of relatively large disks (of tens 
of A U or more) only. Our i nitial rotation is relatively 
fast: Goodman et al. (1993) estimated observationally 
that the ratio ot rotational to gravitational energies for 
dense NH3 cores is typically 0.02, which would imply 
~ 3 X 10'^ cms~^ for our setup. Our choice of a faster 
rotational speed is conservative in that it is harder to 
remove a larger initial angular momentum through mag- 
netic braking (see also Mell on fc Li||2008 ) . 

The calculations arc done with a new version of Zeus, 
dubbed "ZeusTW" . ZcusTW, written in idiomatic For- 



tran95, is base d on the "Zeus36" code ( KrasnopolskyJ 
Li, & Blandford 1999, 2003), itself derivedTrofn Zcus3U 



(L(JA version 3.4.2: |Ciarke, Norman, fc Fiedler||1994| ); 
Zeus36 is parallel by domain decomposition, and utihzes 
dynamic memory for its field and grid arrays. A small 
memory pool of temporary arrays replaces the worker 
arrays of ZeusSD, increasing programming fiexibility at 
essentially zero runtime cost. ZeusTW adds to Zeus36 
the ability to solve many non-ideal MHD problems in an 
explicit form, covering the ohmic. Hall, and ambipolar 
diffusion terms. To treat the ohmic term releva nt for this 
paper, wc used a resistivity algorithm based on Fleming,' 
Stone, fc Hawley^ (2000) , which includes subcycling. Wc 



tested the code by dittusing an initial Gaussian profile in 
Cartesian geometry, in one, two, and three dimensions. 
We also tested slightly different forms of the algorithm 
(such as changing the operator splitting by calculating 
the current density J before or after the magnetic field B 
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Figure 1. Snapshot of the rotationally supported disk formed in 
the absence of magnetic braking at a representative time t = 10^-^ s. 
Logarithm of density (colormap); radial infall sonic transition Vr = 
—a = —2 X 10*cms~^ (white line); poloidal velocity field (arrows, 
scale can be estimated from the radial infall component at the white 
line); and level of Keplerian support \v^/vk\ (dark gray dashes at 
90% and solid line at 100%), where vk = (GA/,/ro)i/2 

is updated), and different subcycling prescriptions (such 
as varying the maximum number of ohmic subcycles from 
1 to 100). The code passed all of these tests. 

2.1. Disk Formation in the Non-Magnetic Case 

We adopt the standard hydro outfiow boundary condi- 
tions at both the inner and outer (radial) boundaries of 
the computation domain (these boundary conditions im- 
pose continuity of material outfiowing the grid, helping 
to reduce artificial refiections of waves at boundaries). 
We use 400 non-uniformly spaced grid points in the ra- 
dial direction, corresponding to a smallest radial grid size 
drmin = 3 X 10^^ cm (or 0.2 AU) and a ratio of adjacent 
zone sizes of 1.017. We use 180 uniform zones in the 6 
direction, corresponding to a minimum polar grid size of 
2.6 X 10^^ cm, comparable to the smallest radial grid size. 
With this setup, we obtain a rotationally supported disk 
of radius w 400 AU and mass rrid ~ O.OIM0 at a time 
t = 10^^ s in the absence of magnetic braking; the time 
is comparable to the duration o f the deeply e mbedded 
Class phase (lAndre et al.||2000l see, however, |Evans et 



al. 2009| who concluded that the duration is longer). A 



snapshot of the disk is shown in Figure [T] The hydro- 
dynamic disk is resolved in the radial direction by more 
than 200 zones. We have experimented with coarser grids 
(e.g., a total of 200 radial zones and 90 angular zones), 
finding qualitatively similar results, which indicates that 
the resolution employed is sufficient. An example of the 
lower resolution runs is shown in ij4] below. 

The disk size at t = 10^^ s can be understood roughly 
as follows. We first estimate the radius of the infall re- 
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Figure 2. Same as Figure mbut for tlie ideal MHD model. The 

field lines are drawn as thin solid lines (20 contour lines of magnetic 
flux, with flux levels quadratically spaced from zero at the axis up 
to the flux value at the outer equator point of each magnetized 
figure). The field component is 470 /^G at the center of the 
figure. No rotationally supported structure forms in this simula- 
tion, although the flow structure is strongly affected by numerical 
reconnection of field lines. 

gion r at time t, and the specific angular momentum £ 
of the material at that radius. Assuming that pressure 
gradients are negligible compared to the central gravity, 
we have the infall speed at any given radius to be near 
the free fall value v{r) = {2GM^/rY/'^. The infaU radius 
is then roughly t ^ r/v ^ r^' V(2G'M,)i/2, which yields 
r ~ 5 X 10^^ cm at t ~ 10^^ s, much larger than zJc (the 
characteristic radius for rotation initial profile change). 
Therefore, the equatorial material at radius r has a spe- 
cific angular momentum ot£ = rv = rx2x 10*cms~^. 
Conservation of angular momentum around a star of 
mass O.5M0 tells us that the specific angular momen- 
tum corresponds to a centrifugal radius of ^ 1000 AU, 
which is about a factor of 2 larger than our numerically 
obtained disk radius ^ 400 AU. We consider the agree- 
ment satisfactory, given the crudeness of the infall esti- 
mate, and the fact that most material on the surface of 
the infall sphere at radius r has a specific angular mo- 
mentum £ less than the above estimate (valid only for the 
equatorial region) , and the centrifugal radius is sensitive 
to £ (as £^). 

2.2. Suppression of Disk Formation in the Ideal MHD 

Limit 

As mentioned in §1, previous studies have shown that 
a moderate magnetic field corresponding to a dimension- 
less mass-to-flux ratio of a few can completely suppress 
the formation of rotationally supported disks in the ideal 
MHD limit (at leas t when the rot ation and magnetic axes 
are aligned; see |Hennebelle &: Ciardi,2009) . Even though 



the mass-to-fiux ratio is ill-defined in our current model 
setup that does not include the self-gravity of the en- 
velope gas, we find that the disk formation can still be 
suppressed if the magnetic field is strong enough. An 
example is shown in Figure [2] where the initially uni- 
form mass distribution in the computational domain of 
the hydro model shown in Figure [T] is threaded with a 
uniform magnetic field Bq = V47r x 10 /iG « 35 /iG 
along the rotatio n axis. The field s t rengt h is in the 
range inferred by Troland & Crutcher ( 2008 ) for a sam- 
ple of dark cloud cores, after correcting tor projection ef- 
fects. For comparison, we also considered a weaker field 
of Bq = Vin X 3 /iG ~ 10.6 fiG and obtained qualita- 
tively similar results. For MHD calculations, we adopt 
the standard outflow boundary condition for the mag- 
netic field at the outer radial boundary, and a torque-free 
(i.e., B^ = 0) outflow boundary conditio n at the inner 
radial boundary, as in Mellon & Li ( 2008 ) . 

The flow structure is strongly attected by numerical re- 
connection of the field lines. The reconnection is an un- 
avoidable consequence of the dragging of the field lines 
into a highly pinched, split-monopole type configu ration 
by gravitat ional collapse in the ideal MHD limit (Galli 
et al.|2006 ). The build up of magnetic flux in the central 



region as a result of continuous mass accretion forces the 
oppositely directed field lines above and below the equa- 
tor closer and closer together, eventually triggering (nu- 
merical) reconnection. It is present in other ide al MHD 
calculations of magnetized core collapse, such as |Mellon| 



fc Li| ( 2008 1 . Because of the (unavoidable) numerical ar- 
tifacts, it IS difficult to quantify the exact strength of the 
magnetic braking. Nevertheless, it is clear from Figure 
[2] that a coherent disk does not exist; it is replaced by a 
set of dense blobs, which are not rotationally supported 
(and may break up asymmetrically in 3D). The promi- 
nent dense blobs and sheets at radii 400 AU and smaller 
are all rotating well below their local Keplerian speed 
(especially the innermost ones), as can be seen from the 
thick dashed and solid dark gray lines at radii greater 
than vj ^ 400 AU in Figure which mark the location 
of \v^\ = 0.9 vk and vk, respectively {vk is the local 
Keplerian speed). The blobs and sheets are supported 
to a large extent by the magnetic tension force from the 
pinched field lines. 

3. INABILITY FOR CLASSIC RESISTIVITY TO ENABLE 
DISK FORMATION 

We now address the question of whether the illustra- 
tive classic resistivity given in Equation ^ can weaken 
the magnetic braking enough to enable the formation of 
rotationally supported disks or not. For this purpose, 
we have implemented into the MHD code ZeusTW an 



Fleming, Stone, & Hawley (2000 


I and Miller & Stone 


( 1997 see also 


Fendt & Cemeljic 


2002 


). In these treat- 



ments, the resistive term of the induction equation is 
included through operator splitting, respecting the con- 
strained transport condition that keeps the field diver- 
gence null to machine round-off. Numerical stability of 
this term limits the timestep At^ to a value oc (Aa;)^, a 
quadratic requirement that can be much more stringent 
than the time AiiMHD oc (Aa;)-^ required for treating the 
hydrodynamics and the ideal MHD terms. Fortunately, 
the resistive term is computationally inexpensive com- 
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Table 1 

Model Characteristics 



Model Bo(pG) ri Disk Notes 



HD 







no 



yes Hydro simulation 



IMHD 

CR 

ER18 

ER18.5 

ER19 

ER19.5 

ER20 



35.4 
35.4 
35.4 
35.4 
35.4 
35.4 
35.4 



no 

yes 

yes 
yes 
yes 
yes 
yes 



no 
no 
no 
no 
no 
yes 
yes 



Ideal MHD 
Classic resistivity given by Eq. ([sjl 
Enhanced resistivity r} = 10^* 
Enhanced resistivity = 3 X 10^^ 
Enhanced resistivity r} = 10^'' 
Enhanced resistivity rj = 3 X 10^^ 
Enhanced resistivity r} = 10^*' 



IMHDw 10.6 
ERlSw 10.6 
ER18.5W 10.6 
ER19w 10.6 
ER19.5W 10.6 
ER20w 10.6 
Note — "Disk" 
disk" at time t - 
cgs units ( cm^ s 



no 
yes 
yes 
yes 
yes 
yes 



no 
no 
yes 
yes 
yes 
yes 



Weaker field, 
Weaker field, 
Weaker field. 
Weaker field. 
Weaker field. 
Weaker field, 



ideal MHD 
: IQiS 

: 3 X 10 
: 10" 
: 3 X 10" 



18 



Tj = 10 



20 



in column 4 refers to "rotationally supported 
= 10^^ s. Enhanced resistivities r] are given in 



pared to the rest of the induction equation, and so it can 
be efficiently subcycled. We have checked that subcy- 
chng speeds up the calculations typically by a substan- 
tial factor and does not change the results significantly. 
We have also tried small changes in the operator split- 
ting scheme, such as altering the order of the ideal MHD 
and the resistivity terms, or adding both terms together. 
These changes had no substantial effect on the results, 
and therefore we chose the scheme that was most con- 
venient for e fficient subcycling, which is one id entical to 
that used in [Fleming, Stone, fc Hawley (2000). 

Figure [3] displays a snapshot of the simulation at the 
representative time 10^^ s. A dense flattened structure 
develops in the inner equatorial region. It is, however, 
fragmented as in the ideal MHD case, indicating that the 
classic resistivity is too small to completely suppress the 
numerical reconnection of field lines. The dense blob at 
radius ^ 400 AU, in particular, is similar to the outer- 
most blob in Figure[2] and is due to an early reconnection 
event. The artificialreconnection events make it difficult 
to carry out detailed analysis. Nevertheless, it is easy 
to show that the dynamics of the inner equatorial re- 
gion is far from that of a rotationally supported disk. 
First, the equatorial material inside ^ 400 AU is often 
counter-rotating, with a rotation speed that sometimes 
decreases toward the center (rather than increasing as 
would be expected for a Keplerian disk: see Figure |4]) . 
Second, the radial motion in the region is mostly super- 
sonic (> 2 X 10^ cms^^). The rapid infall is an indication 
that the dense, flat tened structure is a magnetically pro- 
duced pseudodisk ( Galli &: Shu||1993 ), rather than a ro- 
tationally supported Keplerian disk. There is one region 
of supersonic outflow at ~ 50 AU. This region, and other 
outflow events that appeared during the simulation, are 
related to episodic reconnection events, as in the ideal 
MHD model. In any case, classic resistivity does not 
appear capable of weakening magnetic braking enough 
to enable the formation of a rotationally supported disk. 
Enhanced resistivity is needed if ohmic dissipation is to 
save the disk. 

4. ENHANCED RESISTIVITY AND DISK FORMATION 

It is well known that anomalous resistivity orders of 
magnitude larger than the classic value is required to 
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Figure 3. Same as Figures |T] and [2] but for the MHD model with 
the classic resistivity given m Equation (jsj. A dense, flattened 
structure exists near the equator, but it is a highly dynamic mag- 
netically produced pseudodisk rather than a rotationally supported 
disk. The flow structure is affected by numerical reconnection of 
field lines. 
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Figure 4. Equatorial infall Vr (solid line) and rotation 
(dashed) for the MHD model with classic resistivity, showing the 
lack of a rotationally supported disk. Note the counter-rotation 
inside 400 AU, with rotation speed often decreasing (rather than 
increasing) toward the origin. The spikes of positive radial velocity 
at 50 are related to episodic reconnection events. 



explain the Il-year sunspot cycle (Parker 2007 §9.3). 
Large anomalous resistivity may also exist m star- 
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Figure 5. Same as Figures [T] and [2] but for Model ER18 with an 
enhanced resistivity rj = 10^*cm^ s" . 
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Figure 7. Same as Figures [T] and [2] but for Model ER20 with an 
enhanced resistivity rj = lO^^cm^ s" . 
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Figure 6. Same as FiguresJT] and [2] but for Model ER19 with ; 
enhanced resistivity r? = lO^^m^ s . 



forming molecular cloud cores at high densities ( Norman 
&: Heyvaerts||1985 ), although its magnitude is uncertain. 
Here we consider the simplest case of a spatially uni- 
form resistivity of a range of values and for two different 
magnetic strength Bq = 35.4 and 10.6 /LtG. The models 
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Figure 8. Same as Figure |7| with reduced space resolution. 



are summarized, together with the ones discussed ear- 
lier, in Table [Tl We will concentrate on three representa- 
tive, stronger held (Bq = 35.4 /zG) models with r/ = 10^* 
(Model ER18), lO^^ (ER19), and lO^Ocm^s"! (ER20), 
respectively. The weaker field models yield qualitatively 
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Figure 9. Equatorial infall speed for Model ER18 {rj = 
10^* cm2 s-l, solid line), ER19 {rj = lOi^cm^s-l, dashed), and 
ER20 (»7 = lO^'^cm^s"^, dash-dotted line) at the representative 
time t = 10^^ s. The free-fall speed is plotted as the lower dotted 
line for comparison. 

similar results, and will be discussed briefly toward the 
end of the section. 

Figures [5}|7] show snapshots of the three models with 
enhanced resistivity, and Figures [9f |ll| display their in- 
fall and rotation speeds and vertical field strength on 
the equator. Additionally, Figure [8] shows the effect on 
Model ER20 of reducing numerical resolution by a factor 
of 2 in both r and 0. Compared to the models without 
resistivity (Model HD shown in Figure IlJ and with the 
classic resistivity (Model CR, Figure |3| , the mass distri- 
bution and velocity field are much smoother, indicating 
that the enhanced resistivity in these models has largely 
eliminated the numerical reconnection. The suppression 
facilitates quantitative analysis of the simulations. 

The mass distributions in the two lower resistivity 
models (Models ER18 and ER19) appear similar. Both 
have a dense equatorial region sandwiched from above 
and below by a more diffuse infalling envelope. The infall 
in the inner part of the envelope (within a few hundred 
AUs) appears to be guided by the field lines. It strikes 
the surfaces of the dense equatorial structure at a super- 
sonic speed, forming two "accretion shocks." These are 
reminiscent of the classic accretion shocks around Ke- 
plerian disks formed in n on-magnetic simulations (e.g., 
Bodenheimer et al. 1990 see Figure [l]). Their produc- 
tion is due, however, mainly to the magnetic field rather 
than rotation. Indeed, the same shocks form even in 
the absence of any rotation. For this reason, we will 
term these shocks "the magnetically induced accretion 
shocks." They surround a m agnetically induced "pseu- 
dodisk" ( [Galli fc Shu| 1993| rather than a rotationally 
supported Keplerian disk. 

Figure |9] shows that the dense equatorial structure in 
Models ER18 and ER19 contains a pseudodisk rather 
than a rotationally supported disk. In both cases, the 
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Figure 10. Same as Figure[9]but for the equatorial rotation speed. 
The Keplerian speed is plotted as the upper dotted line for com- 
parison. 
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Figure 11. Same as Figure |9] but for the vertical magnetic field 
component Bz = —Bg on the equator. 



infall is supersonic (> 2 X lO'^cms"!) within lO^AU of 
the central object, although there is significant difference 
in the infall speed within ~ 200 AU: the equatorial ma- 
terial in the less resistive model ER18 falls inward much 
more slowly than that in the more resistive model ER19. 
The prominent deceleration of infalling material around 
~ 200 AU in Model ER18 is reminiscent of the ambipo- 



cussed in 


Li & McKee ( 1996 see also iCiolek & Konigl 


19981 Li ; 


[m\ 'Rrasnopolsky & Konigl i>00i>| 'i'assis k 


Mouschovias..2007, and especially Figure 3 of Mellon & 
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Li|2009 ). The physical origin is also similar. Just as am- 
bipoiar diffusion, the ohmic dissipation enables the field 
lines to diffuse relative to the material that plunges into 
the central object. The net effect is the formation of a 
strong (poloidal) field region at small radii, where the 
inward advection of the field lines by infalling material 
is balanced by the resistivity-induced outward magnetic 
diffusion. A relatively small resistivity such as that in 
Model ER18 allows for relatively larger magnetic gradi- 
ents and Lorentz forces, which decelerate the infalling 
material more strongly. When the resistivity is small, 
the magnetic field remains a significant barrier to proto- 
stellar accretion. 

The magnetic barrier weakens as the resistivity in- 
creases. In Model ER19, which is ten times more re- 
sistive than Model ER18, the ohmic dissipation-induced 
deceleration has all but disappeared. This is because 
the higher resistivity can support only a smaller cur- 
rent, making it hard to effectively "trap" the magnetic 
flux that is dragged to the small radii by the collaps- 
ing flow. Nevertheless, the infall speed is considerably 
smaller (by about a factor of two or soVthan the free-fall 
speed, which is also plotted in Figure |9] for comparison. 
The sub-free-fall motion comes about because the equa- 
torial infall is retarded by a combination of centrifugal 
and magnetic forces, although the magnetic forces ap- 
pear to be the more important co ntributor, since the 
rotation speed shown in Figure 10 is about 1/3 of the 
Keplerian speed or less (yielding a centrifugal force an 
order of magnitude weaker than the gravity). 

Note from Figure [TO] for the least resistive model 
ER18, the equatorial rotation speed within 200 AU 
of the origin decreases, rather than increases, toward the 
center. Obviously, the decrease is due to magnetic brak- 
ing, which is apparently so efficient as to produce a region 
of counter-rotation at ~ 100 AU. The counter rotation 
is reminiscent of the case of magnetized core collapse in 
the presence of ambipolar diffusion with a relat ively high 



cos mic ra y ionization rate 10 s ^ studied by Mellon & 



Li ( 2009 see the bottom panel of their Figure 4). In both 



cases, the efficient braking is due to the strong poloidal 
field trapped by infall at the small radii. It removes es- 
sentially all of the angular momentum of the equatorial 
material that crosses the inner boundary of the compu- 
tation domain. In the more resistive model ER19, the 
equatorial rotation speed inside ^ 200 AU is initially 
kept nearly constant, before increasing rapidly. The peak 
rotation speed remains well below the local Keplerian 
speed, indicating that a rotationally supported disk is not 
formed. We conclude that the dense, flattened, equato- 
rial structure produced in both Models ER18 and ER19 
feature dense pseudodisks, but no rotationally supported 
dense disks. These pseudodisks are thinner than the Ke- 
plerian disk shown in Figure [T] for the hydro case because 
of vertical compression by pinched field lines, and smaller 
column densities (due to the fast infall rate inside the 
pscudodisk), corresponding to reduced scale height for 
a given gas density. The absence of a rotationally sup- 
ported disk in these two cases (which have resistivities 
much larger than the classical value) supports our re- 
sults in §2.2 and 3 that disk formation is suppressed by 
magnetic braking in the ideal MHD limit and in the pres- 
ence of the classical resistivity, even though these results 
were affected by numerical reconnection. 



The most resistive model ER20 differs drastically from 
the lower resistivity cases. It produces a much denser 
disk of radius ~ 150 AU at the representative time t = 
10^^ s (see Figure [7|, which turns out to be rotationally 
supported. The strongest evidence for rotational support 
comes from Figure [TOj which shows that the rotation 
speed at the equator is nearly identical to the Keplerian 
speed inside ~ 150 AU. The support is also evident in 
Figure [9] which shows that the infall has nearly stopped 
at smallradii (the spike near the inner edge of the com- 
putation grid is likely related to the outflow boundary 
condition, which is not ideal for a rotationally supported 
disk; it was present in the hydro case as well). The 
thick dashed and solid dark gray lines in Figure [fj which 
mark the location of \v^\ = 0.9 vk and vk, respectively, 
also make the case. The disk has a mass ~ O.OIM0 at 
t — 10^^ s, corresponding to an average accretion rate 
(from the envelope to the disk) ~ 5 x 10~^Mq yr~^. As 
in the pure hydro case (§2.1), little of the disk material 
passes through the disk onto the central object in our 
axisymmetric simulations. Gravitational torques may be 
needed to drive the further evolution of such disks (in 
three dimensions). 

Another feature of Model ER20 is that the poloidal 
field lines pass through the equatorial region with only 
a moderate bending. It is of course a consequence of 
the large resistivity, which allows the field lines to diffuse 
outward easily relative to the infalling material. This is 
especially true inside the inner, Keplerian disk, where the 
combination of slow radial motion and large resistivity 
means that the magnetic field must be nearly vertical 
and uniform (see the dot-dashed line in Figure 11). A 
consequence of the nearly uniform field is that there is 
little vertical compression by the magnetic field, yielding 
a Keplerian disk that is thicker than the magnetically 
compressed pseudodisks in the less resistive models (see 
Figures [5] and |6 1 . 

Outside the I<^eplerian disk lies a supersonically in- 
falling pseudodisk. It is linked to the surrounding en- 
velope by magnetic field lines; they are separated by the 
magnetically induced accretion shocks. Upstream of the 
shock, the low density envelope collapses more or less 
along the field lines, indicating that flux- freezing is ap- 
proximately valid, because the field lines there are rela- 
tively straight and there is little current to be dissipated 
by the (large) resistivity. Downstream of the shock, the 
dense pseudodisk collapses horizontally, nearly perpen- 
dicular to the (poloidal) field lines. At the relatively late 
time shown in Figure n\ the field lines inside the pseu- 
dodisk are nearly helafixed in space, with the inward 
advection due to infalling material balanced by the out- 
ward resistivity-driven diffusion. 

The role of resistivity i n e nabling disk formation can 
be seen clearly in Figure |11[ where the vertical compo- 
nent of the magnetic field on the equator, B^, is plotted 
as a function of time for all three models with different 
enhanced resistivities. The equatorial infall in the pseu- 
dodisk tends to compress the poloidal field lines, leading 
to a sharp increase in the poloidal field strength at small 
radii. Resistivity, on the other hand, tends to smooth out 
the field distribution. The net effect is that the poloidal 
field strength at small radii depends strongly on the mag- 
nitude of the resistivity, being larger for smaller resistiv- 
ities. Since most of the magnetic braking occurs at rela- 
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tively small radii, where the infalling material attempts 
to spin up the fastest to conserve angular momentum, 
the stronger poloidal fields in the two less resistive cases 
(models ER18 and ER19) produce a stronger braking, 
which is enough to suppress the formation of a rota- 
tionally supported disk. The most resistive case (model 
ER20) has the weakest poloidal field at small radii (where 
magnetic braking matters most in the model) which, to- 
gether with a large slippage between the field lines and 
matter in the azimuthal direction (which makes it hard 
to generate a toroidal field) , makes the magnetic braking 
ineffective and disk formation possible. That the disk 
in model ER20 forms as a result of enhanced resistiv- 
ity rather than numerical diffusion is supported by the 
fact that a nearly identical disk is formed in the lower 
resolution run (see Figure |8| . 

One can roughly estimate the size of the resistivity- 
enabled disk as follows. In order for the formed rotation- 
ally supported disk to survive against magnetic braking, 
the disk rotation must not twist up the poloidal mag- 
netic field that threads the disk too strongly. This is sat- 
isfied when the magnetic diffusion speed becomes com- 
parable to the Keplerian speed, which leads to the con- 
dition rj ~ VfcH, where H is the vertical scale length 
over which the toroidal magnetic field varies. If we take 
H to be the disk pressure scale height (ignoring possible 
magnetic compression of the disk) , and use the standard 
result for a thin-disk H/r ~ a/vK (where r is the ra- 
dius and a is the isothermal sound speed), then rj ra 
or r ^ 7]/ a ^ 3.3 x 10^ AU (for r] = lO^^cm^s"^ and 
a — 2x lO^'cms^^). This is about a factor of two larger 
than the disk radius that we found numerically in Figure 
[71 We should caution the reader that the above estimate 
did not take into account the detailed magnetic geom- 
etry and braking efficiency, which depends on the field 
strength (see below) . For a weak enough field, a disk can 
form independent of the magnitude of the resistivity. 

Besides the three representative models discussed 
above, we have carried out two additional simulations 
with r; = 3 X lO^^ (Model ER18.5) and 3 x 10^^ cm^ s''^ 
(Model ER19.5). As one would expect based on the 
results for Models ER18 and ER19, no rotationally 
supported disk forms in Model ER18.5. There is a 
small (~ 40 AU in radius) rotationally supported disk 
in the more resistive Model ER19.5 at the fiducial time 
t — 10^^ s. The disk shrinks with time, however. By 
t = 2 X 10^^ s, its radius decreases to ^ 25 AU; it may 
disappear altogether at a later time. The trend indicates 
that 77 = 3 X 10^^ cm^ s~^ is probably close to the critical 
resistivity 77c needed for disk formation. 

The value for the critical resistivity rjc depends on the 
initial field strength Bq. This is to be expected since, in 
the limit of infinitely weak field, rotationally supported 
disks can form without any (enhanced) resistivity at all 
(i.e., 77c — 0). For a moderately magnetized dense core 
of Bq = 10.6 /zG, we find that the transition between 
the formation of a rotationally support disk and its sup- 
pression occurs between rj = 10^^ (Model ER18w) and 
3 X lOi^cm^s-i (Model ER18.5w). This value of Bq is 
probably close to the lower limit to the field strength 
in dense cores, judging from the fact that the median 
field strength for the more diffus e, cold n eutral medium 
(CNM) of atomic gas is - 6 ^iG ( jHeiles fc'Troland||2005| 



and that the directly measured line-of-sight component 
of the magnetic fiel ds in a number of dense co res is of 
this order or higher (Crutcher & Troland 20081); the full 



strength of the core magnetic field is likely significantly 
higher. It is therefore reasonable to expect the critical 
value rjc for disk formation in cloud cores magnetized to 
a realistic level to lie somewhat between 77c ~ 3 x 10^^ 
and ~ 3 X 10^^ cm^ s""'^. In what follows, we will take 
rjc — lO^^cm^s^^ as the characteristic critical value, 
with the understanding that it depends somewhat on the 
field strength, and can be uncertain by a factor of ~ 3 in 
either direction. 

5. DISCUSSION AND CONCLUSION 

Our most important result is that classical resistivity 
alone appears unable to weaken magnetic braking enough 
to allow a rotationally supported disk to form for a re- 
alistic level of core magnetization, and that an enhanced 
resistivity of order 77c — 10^^ cm^ s~^ is needed. 

A large, enhance d resistivi ty (of order 2 x 10^° cm^ s~^) 
was advocated by |Shu et al. (2006) to solve the so-called 
"magnetic flux problem" m star formation. They intro- 
duced the concept of the Ohm sphere, with a radius: 



roh, 



2GAL 



ill) 



where Af* is the central stellar mass. It can be obtained 
by setting the characteristic magnetic diffusion speed 
= rj/r equal to the free-fall speed vg — (2GM^,lr)^l'^ . 
In the limit of a weak magnetic field (so that the mag- 
netic forces d o not affect the a ccretion dynamics) and 
zero rotation, Shu et al. (2006) showed that the mag- 
netic field becomes more or less uniform inside the Ohm 
sphere; the sphere may be thus regarded as the region 
of complete magnetic decoupling. If a rotationally sup- 
ported disk were to form inside the Ohm sphere, it would 
be protected against strong braking by the decoupling. 

If the Ohm sphere must be larger than the rotation- 
ally supported disk, then to form a disk of 10^ AU in 
radius aro und a ha lf-solar mass star, the resistivity must 
be 770 — \J2GM^r^ ~ 4.5 x 10^*^ cm^ s~^. Our actual, nu- 
merically obtained, value of resistivity for disk formation 
is one to two orders of magnitude smaller, depending on 
the magnetic field strength. The discrepancy indicates 
that the simple considerations leading to the Ohm sphere 
need to be modified. One such modification is that the 
(radial) magnetic diffusion speed is expected to be 
larger than ?7/r, because the field lines bend sharply in- 
side the pseudodisk, with a field variation scale A much 
less than r (see Figs. [5] and |6]) . The other is that the 
infall speed can be significantly below the free-fall value 
(see Fig. [9]). Both modifications make it easier for the 
field lines to diffuse out, lowering the critical value for 
disk formation compared with the simple estimate based 
on the concept of Ohm sphere. 

Nevertheless, the critical resistivity for disk formation 
is still a few orders of magnitude higher than the classi- 
cal value. H ow such an en hanced resistivity may arise is 



unclear (see [Parker 1979). One possibility is turbulent 
resistivity, which has been discussed, for exa mple, in the 
context of magnetized accretion disks (e.g., Lubow, Pa- 
paloizou, fc Pri ngle'1994{ [Guan fc Gammie|2009[ ) . There 
are at least two possible origins tor turbulence on the disk 
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formation scale (of 10^-10'^ AU). First, the dense cores 
that collapse into stars are observed to have non-thermal 
motions on the tens of thousands of AU scale, which can 
be interprete d as a subsonic turbule nce of order half the 
sound speed (Bcrgin & Tafalla 2007) . It is likely that the 
turbulent motions would continue down to smaller scales 
as the core collapses. Second, the collapsing flow can be 
significantly retarded by magnetic forces (see, e.g., the 
infall speed for Model ER18 in Figure |9]). The magnetic 
support of the collapsing material against gravity may 
lead to interchange instability in 3D, which could de- 
velop into a turbulence, and might itself lead to an effec- 
tive turbulent resistivity. Whether turbulent resistivity 
can enable the formation of rotationally supported disks 
or not requires high-resolution 3D s imulations, perhaps 
along the line of Kowal et al. (2009|). Another possibil- 
ity is anomalous resistivity associated with, for example, 
current-driven instabilities (Norman fc Hcyvaerts, 1985( ). 
Whether these and other ennancements of resistivity can 
enable disk formation remains to be determined. If not, 
other means of enabling disk formation, such as the Hall 
effect and outflow stripping of the protostellar envelope, 
need to be explored. 

Although our idealized calculations indicate that a con- 
siderable enhancement in resistivity is needed if ohmic 
dissipation alone is to enable the formation of rotation- 
ally supported disks during the protostellar accretion 
phase of low-mass star formation, the exact value for the 
required enhancement is uncertain, and can be better de- 
termined with future refinements. One refinement is to 
include the self-gravity, which would also allow us to fol- 
low the formation and evolution of the dense magnetized 
core to the point of point-mass formation. The prestellar 
core evolution would provide a more self-consistent ini- 
tial condition for the disk formation calculation during 
the protostellar accretion phase. An integrated treat- 
ment of both core formation and protostellar accretion 
may enable us to determine the distribution of the matter 
surrounding the dense core, which can aid the core mate- 
rial in removing the angular momentum from the forming 
disk, especiall y if the core is hig hly flattened along the 
field hnes (|Basu fc Mouschovias||1994, Krasnopolsky fc] 



|Konigl|2002[ see e.g., |'l'assis et al.||Wgl tor statistical ev- 
idence of core flattening). Another refinement, which we 
plan to do, is to study the combined effect of resistivity 
and ambipolar diffusion on disk formation. In addition, 
it would be desirable to find a better way to quantify the 
numerical reconnection, originated physically from the 
strong collapse-induced magnetic pinching, that affects 
the ideal MHD and low resistivity simulations. We con- 
clude that how protostellar disks form around accreting 
protostars remains a mystery. 
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